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ABSTRACT 


The presence of soot on the fuel side of a diffusion flame results in significant radiative 
heat losses. The influence of a fuel side heat loss zone on a pure diffusion flame established 
between a fuel and an oxidizer wail is investigated by assuming a hypothetical seek 2 heat loss 
profile. The intensity and width of the loss zone are parametrically varied. The loss zone is 
placed at different distances from the Burke-Schumann flame location. The migration of the 
temperature and reactivity peaks are examined for a variety of situations. For certain cases the 
reaction zone breaks through the loss zone and relocates itself on the fuel side of the loss zone. 
In all cases the temperature and reactivity peaks move toward the fuel side with increased heat 
losses. The flame structure reveals that the primary balance for the energy equation is between 
the reaction term and the diffusion term. Extinction plots are generated for a variety of 
situations. The heat transfer from the flame to the walls and the radiative fraction is also 
investigated, and an analytical correlation formula, derived in a previous study, is shown to 
produce excellent predictions of our numerical results when an 0(1) numerical multiplicative 
constant is employed. 


2 



1 . 


INTRODUCTION 


The interaction between the structure of a diffusion flame (DF) and the flame radiation 
is quite complex. Soot is formed and oxidized in a diffusion flame as a consequence of a variety 
of physical and chemical processes. There are considerable uncertainties in the description of 
soot processes in a flame and the soot evolution mechanisms are not completely understood. 
Hence, the solution of the complex problem of diffusion flame - soot radiation interaction is very 
involved. The energy, species and soot volume fraction equations are all coupled and contain 
nonlinear source terms. In the present work we investigated the influence of a simple and 
contrived heat loss profile on a pure diffusion flame established between two diffusing walls of 
fuel and oxidizer. 

The influence of heat transfer by radiation on flames has received significant attention 
in recent years [1-7]. Thermal radiation from a flame can occur from (1) radiation from the 
combustion gases at high temperature and (2) radiation from combustion generated particulates, 
i.e. , soot. According to the calculations of Grosshandler and Modak [8] for soot volume 
fractions greater than 10' 7 soot radiation is dominant. A review of pure diffusion flames without 
heat losses is presented first. In the following sections we define the problem geometry, 
describe the particular form of the heat loss profile used, formulate the conservation equations, 
briefly indicate the numerical method used and discuss the results. 

In an ordinary diffusion flame the characteristic flow time is much greater than the 
characteristic chemical reaction time. This implies that the chemical reaction is much faster than 
the transport of species to the flame unless the flame is near or approaching the extinction stage. 
A pure diffusion flame is established when both oxidizer and fuel are transported to the flame 
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by means of diffusion only. No convective flow is present. Some important characteristics of 
pure diffusion flames have been discussed in [9], including an analysis of the detailed nature of 
the temperature and reaction rate profiles. It was observed that the maximum of the reaction 
rate profile usually will not coincide with the temperature profile maxim. The only exception 
is the symmetric flame for which the overall stoichiometric coefficient, <K = v Y^Y^, equals 
unity. This study shows, for a fuel-rich flame, that i.e. , the peak of the reaction 

rate profile (Zp lies between the Burke-Schumann flame location (Zp and the peak of the 
temperature profile (Z t ) for fuel rich conditions. For oxidizer-rich conditions, Z.<.Z r <>l f 
It may be argued that in thin-flame limit all diffusion flames are "pure" diffusion flames because 
the mixture fraction transformation discussed in Williams, chap. 3 [10], produces an equation 
resembling T^-lVZl' 2 ^, where [ VZ [ is the magnitude of the mixture fraction gradient 
perpendicular to the flame. However, | VZ | depends strongly on the heat and fluid flow 

conditions and in effect introduces a new parameter that must be accounted for in a complete 
analysis. Hence, though the value of |VZ,| (i.e., |VZ| evaluated at the flame sheet) may 
be buried into a suitably redefined Damkohler number, it must of course be resurrected when 
later conducting a full examination of the problem. 

Our goal in this work is to describe the response of a diffusion flame, when there is a 
zone of radiant soot-generated energy losses nearby, through the examination of a simplified 
model. A previously-generated theoretical correlation will be tested, and we shall, in addition, 
attempt to produce practical correlations of the total heat flux from the flame, the total 
(conductive plus radiative) energy flux to the surface, the drop in flame temperature due to 
radiant loss, etc. We shall not develop the full correlations here, but we indicate their dominant 
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behaviors. The complete correlations require numerical examination of all possible cases 
something which was beyond the scope of this study. 


2. THE MODEL 

Figure 1 schematically depicts the problem geometry. The fuel wall and the oxidizer 
wall are located at ;c = 0 and x=L , respectively. Both walls issue diffusive fluxes of the 
respective constituents. A diffusion flame is established between the two walls. A soot layer 
is assumed to exist on the fuel-side of the DF, consistent with experimental observations [11]. 
The walls have the ambient temperature T 0 . There is no fuel on the oxidizer wall and no 
oxidizer on the fuel wall, the only possible other species at the walls being an inert element. 
The fuel and oxidizer mass fractions at the respective walls are specified to have values Yfr and 
Yoo as shown in Figure 1 . 

The combustion reaction is assumed to be a global, one-step chemical reaction of the 
form F+vO-(l +v)P, where F denotes the fuel and 0 denotes the oxidizer. Methane is 
nominally the fuel under consideration and oxygen is the oxidizing specie, although real 
methane-oxygen reactions require of the order of 100 reaction steps and individual property 
choices for the separate species. The stoichiometric fuel-oxidizer mass ratio, v, is four for the 
methane-oxygen combustion reaction. A suitable set of parameter values must be used to 
generate a reasonable range of Damkohler number and flame temperature values. The adiabatic 
flame temperature is given by Tf=T 0 +Q F Y FF l[C p (\+§)], where Q f is the heat release per 
unit mass of fuel from the combustion reaction and C p is the specific heat of the mixture. 
The global stoichiometric coefficient, denoted by 4>, is given by vY FF }Y ^ However, the 
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use of the above formula produces unrealistically high adiabatic flame temperatures. Thus, the 
above equation for T f was modified to produce a practical and reasonable range of adiabatic 
flame temperatures. 

A set of realistic hydrocarbon combustion flame temperatures was obtained from the 
work of Wichman [12] for the analysis of flame spread over thermoplastics. The idea there was 
that the fuel mass fraction cannot reasonably be determined at the surface but a more-or-less 
generic flame temperature can still be calculated. This flame temperature varies only with the 
free-stream oxidizer mass fraction Y^. The value of Y fr , i.e., the fuel mass fraction in 
the hypothetical fuel stream for our present calculations was assumed to be 0.85. The oxidizer 
mass fractions ( Y ^ and the flame temperatures ( 7 are tabulated in Table 1. 

A fourth order polynomial was fitted to the above data to obtain 

7^486.66 + 12230.857^-25728. 647^+25360.027^,-9323.07^. (1) 

Using this expression, we generate points for a (7^7^.) plot. 

We calculate Q F by using the relation Q F ={J f -T^)C p {\ t§)/Y ff for 
7^=0.211, y^=0.85 and T g =29BK. The calculated value of Q F is 1 1959.43 kJ/kgK. 

We now introduce a modified formula for the calculation of the flame temperature, 


„ ^ . QMYoJ 

f ° C,( i + 4» 


(2) 


Next we calculate the values of the modification factor, A*oo)> by us ' n S above expression. 
The calculated value of Q F and the (7^,7,) data obtained using equation (1) were utilized 
for this purpose. The A^oo) data were as follows, 
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f-\.25t^{-2.99Y^0.33 (3) 

Finally, we use the above expression for fiY^ to calculate for any set of Y ^ and Y FF 
values in equation (2). A plot of T f versus Y ^ is shown in Figure 2. The Y FF values 
corresponding to the different curves in the plot range from 0.25 to 1.0. The lowest curve is 
for Y ff = 0.25. The curves above are for Y FF - 0.30, 0.35, 0.40, etc. We note that for 
T^O.25 and 0.30 the peak flame temperature does not occur at 7^ = 1; there is a slight 
local maximum in the range 0<F oo <l. For this reason we shall not use these curves. 
However, for higher values of Y FF we do obtain temperature profiles whose maxima occur 
at 7^= 1. These profiles will be used. Also, we note that since we are interested in the effect 
of soot radiation on diffusion flames, we are not concerned with low values of which do 

not produce significant amounts of soot. Hence, in our analysis, values of 0.30 and 

lower are not used with reasonable physical justification. 

The parameter values in the work of Tzeng et al. [13] were used in this article. The 
important values are shown in Table 2. 

Here we write the equations and boundary conditions for conservation of energy and 
species. The energy conservation equation is 

pC e lT*uT,MlTX+Q F Ur~- ( 4 > 

w'ith boundary conditions T\x=0)=T g and T(x=L)=T o , where T g is the temperature at the 
fuel and oxidizer walls, assumed to be 298 K. Here T is the temperature, p is the density, C p 
is the specific heat of the mixture, X is the thermal conductivity and u is the velocity. The 
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volumetric radiative heat loss term is -dqjdx with units Wjm'\ The heat release due to 
combustion is Q F and u F is the reaction rate term. An Arrhenius type expression was 
assumed for the one-step irreversible reaction, so the reaction rate is u F = pAY 0 Y F exp-E/RT. 
The quantity A denotes the pre-exponential factor. The thermal conductivity is X. The 
oxidizer mass fraction equation is 

P[ Y o^uY 0 )=(P d o Y o)x- v ^p ( 5 ) 

with boundary conditions T o (x=0)=0 and Y 0 (x=L)=Y 00 . Here, D Q is the mass diffusivity 
of the oxidizer. Similarly the fuel mass fraction equation is 

p[Y F +uY F MpD F Y F ) r u p ( 6 ) 

I i j 

with boundary conditions Y F (x=0) = Y FF and Y f (x=L)= 0; D f is the mass diffusivity of the 
fuel. 

The above equations are now transformed to a mass coordinate system. The transformed 
coordinate is 2=1-5/^ where s-J^pdx and s 0 =j^ L pdx. We note that Z- 1 when ,r = 0 
and Z=0 when x= 1. The coordinate Z happens to be identical to the mixture fraction coordinate 
for our simple problem. The following expressions hold for the above transformation: 


and 


a 

= 

--((pu) 0 -(pu) s )^((pu) 0 -(pu) L ) 

a 

a 

+ — 

dr 

I 

V X 2 u ^ 

J 0 .. 

dZ 

i * 


( 7 ) 

z 


_a = __p_a 

dx t s 0 d z, 


( 8 ) 


Since we are considering a pure diffusion flame, (pu)_ r _ 0 =(pu) i _, =0, i.e., there are no 
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convective flows from the walls. Application of these operators to the energy and species 
equations and assuming p A, p 2 D 0 , p 2 D F are constant, assuming the Lewis number is 
unity and D Q =D F and then introducing the nondimensional variables t= (7- 7^/(7^,- 7^, 
y 0 ~ Y J Y oo’ y f~ Y f^ Y fp , s~^/(Pc/') gives 


1 - N R dq 

'7-i'zz +<}&+■=* * 


-2 

S n 


s 0 dz 


(9) 


y^- 2 yo n -^ r > 


(10) 



(ii) 


where s 0 =sJ{p g L)= J l pdx where p=p/p 0 and x=xjL. Also, r is the nondimensional 
reaction rate and N R is a radiation number [14] evaluated as the ratio of the reference 
radiative and conductive fluxes, given by N R =q iO J(\ 0 (T / -T o )IL). The thermal conductivity 
at the reference condition is denoted by A 0 . The quantity D is the Damkohler number given 
by where the reference diffusion time is t r ^=L 1 ja Q and the characteristic chemical 

time is t ekem = ll[AY 00 txp(-E/RT^]. The nondimensional quantity q R is given by . 

where , is a reference radiative heat flux. The nondimensional heat release, Q p is 
given by Q F Y Ff j[C i> {J^-T^] and equals (l+4>) since the adiabatic flame temperature is 
defined as Tf=T o +Q F Y FF J[C (1 +4>)]. We note that in the prefactor multiplying the reaction term 
of equation (9) we do not utilize the temperature correction discussed previously. In addition, 
we have defined t=t(t. The nondimensional reaction term, r, is written in the form 
r= )'o>’f ex P[ _ P(l -x )/(l -a (l -T ))]> where a=l-TJT / and p=£a/(£ u 7p; £ is the activation 
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energy and R u is the universal gas constant. The quantity P is known as the Zeldovich 
number. 

2.1 Infinite Reaction Rate (DRJR) Solutions as Initial Conditions 

Equations (9)-( 11) are the governing conservation equations for t, y Q and y F for 
the case of finite rate chemistry. The equations become much simpler when the reaction rate 
is infinite. In this case, all fuel reaching the flame surface is consumed instantaneously, and 
similarly for the oxidizer. Thus no fuel exists on the oxidizer side and no oxidizer exists on 
the fuel side, i.e., y Q y f .= 0 on both sides of the flame. The energy equation can now be 
solved in two adjacent domains, the oxidizer side (OsZsZp and the fuel side (Z^sZs 1) of 
the flame. The flame location is designated by Zj. For our simple problem the mass 
coordinate Z coincides with the mixture fraction, a conserved scalar: this is a rare special case. 
In the absence of radiative losses the steady state energy equation becomes t^O. Since, at 
the flame t = 1, the solution of the steady state energy equation ("^=0, since N^=0 and 
the reaction term can be excluded) for the infinite reaction rate (1RR) situation gives t=Z/Z / 
for 0 <.Z<,Zj and t=(l -Z)/(1-Zp for Z^sZsl. Similarly, y 0 and y F can be solved 
for the infinite reaction rate situation and we get y 0 =(l -Z)-(l -Z^x for 0 <,Z$Z f and 
y F =Z-Zp. for Z^sZsl. 

Next, we must evaluate Z. the coordinate location of the IRR flame. As mentioned, 
Z is the mixture fraction coordinate, defined as Z=(4>y r + l-y 0 )/(4) + l). At the flame, y Q and y F 
are both zero so that 2p= 1/(4) + 1). With knowledge of Z p the nondimensional temperature 
and species equations can all be determine exactly. The profiles so obtained are used as initial 
profiles for the numerical solution of the transient conservation equations (9)-(ll). 
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2.2 Simple Heat Loss Profile 

As shown in Wichman and Ray [14] the simplest model heat loss profile is the "top-hat" 
profile used therein. Because of the discontinuous derivatives at the edges of the top-hat profile, 
it is not as convenient for numerical reasons as a smooth and continuous heat loss profile. For 
primarily this reason, the profile that we shall use here is of the form of a seek ' in mixture 
fraction space, viz., 


~ =sech\B{Z-Z R , )). (12) 

dZ 

The location where the maximum of -dq^dZ occurs is denoted by Z R . Figure 3 illustrates 
the nature of variation of the heat loss profile in mixture fraction space. We note that the 
maximum value of the profile is unity. We define the Z locations where the value of the 
function is 1 % of its peak value as the two tails of the function, located respectively at locations 
we presently call Z R _ and Z R- , with Z R _<Z R _. The maximum of the -dq^dZ profile 
occurs at Z R =(Z R _+Z R d/ 2. The width of the loss zone is defined to be A Z R =Z Rt -Z R _. The 
separation distance of the loss zone form the location of the ideal Burke-Schumann flame Z f 
is given by A =Z R _-Z f In the subsequent analysis, we shall vary the thickness A Z^ as well 
as the separation distance A, in order to study the influence of the loss zone on the flame 
structure. The thickness of the loss zone can be chosen by selecting different values of the 
parameter B in equation (12). 

From equation (9) we note that the radiative loss term is given by (NJs^dq^dZ, and 
hence, another important way to modify the loss term is to experiment with the value of its 


11 



amplitude, N^s 0 . We can choose different values of N^, the radiation number. The 
quantity given by ^(p^L) is a by-product of the solution and for this reason is 

evaluated at each time step. 

We recall that for the top-hat profile the integrated heat loss is 

j^(N f Js^(U(Z R _)-U(Z R J)dZ=N J Js 0 &Z Jt . In this case the integrated heat loss is given by 
5 0 l J o N f £ech 2 [B(Z-Z R )]dZ=N f J(Bs l JtaiihB(l -Z R )+ianhBZ R ], For large B this simplifies to 
INjJSo/B+OiB' 2 ), showing that the top-hat loss zone thickness A Z R corresponds to 2 IB, 
or B=2/AZ R . Consequently, in analytical formulae for the top-hat profile (see ref. [14]) we 
can substitute for t\Z R the value 2 IB in order to test their correspondence to the seek ’ profile. 

3. NUMERICAL SOLUTION 

Equations (9), (10) and (1 1) were numerically solved with a the finite difference scheme. 
The nonlinear source terms were linearized using Newton’s method. For each time step, 
iterations were used until the sum of normalized residuals became smaller than 1x10^. The 
transient conservation equations were integrated to steady state. 

We utilize the ideal gas law to derive a relation between the temperature and density of 
the system. We can write pV-(m/W)R u T where p is the density and R is the gas constant for 
the mixture, given by R=RJW. If we assume PqRT q to be the constant pressure of the 
system, then the introduction of a=\-TJT f and x =(T-T ( )l{T f -T i ) results in the following 
important relation: 


P = 


1-g 

l-a(l-t)' 


(13) 
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We observe that when the temperature is that of the ambient, i.e. , T=T g =298K, then t=0 

and p = l, i.e., p = p 0 . At the flame temperature (7jj) the nondimensionai density is 

p = l-a and consequently p=(l-a)p 0 . 

Equations (9)-( 11) indicate that in order to solve the r, y Q and y F equations, we 
need to evaluate s 0 at every time step. The quantity s Q enters the analysis by virtue of the 
coordinate transformation Z=1 -J/J Q . By differentiating both sides of this relation, we obtain 
dZ/dx = -p/s 0 , since s=J*pdx and x=xjL. Using the transformation relation between Z 
and x, subject to the above mentioned boundary conditions, we get 
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and the relation between the x and Z coordinates can be written as 


(14) 


f (1/P )dZ 

x= — - ■ (15) 

f (1 l?)dZ 
j o 

On obtaining the solutions for ~,y 0 and y ^ equation (15) is used to transform the solutions 
back to the physical coordinate x. Thus, J 0 can be evaluated once the p distribution is 
known. The normalized density p can be related to the t distribution by virtue of equation 
(13). Hence, the expression for s 0 can alternatively be written as 


S n = - 


J^dZ+(c/(l -a)) 


(16) 


13 



Hence, on obtaining the t profile we can determine the quantity J Q . When the temperature 
throughout the domain is the same as the ambient temperature then t=0 everywhere; 

by using equation (16), we obtain s 0 = l. If we next assume that the temperature everywhere 
in the domain is the adiabatic flame temperature T p then t = l and s 0 =(l-a). Since the 
minimum and maximum values of temperature are T g and T, respectively, the quantity 
J 0 must obey the limits (l-a)^s 0 sl. 

4. RESULTS AND DISCUSSION 

Figure 4 depicts the nondimensional temperature, t, plotted as a function of the 
mixture fraction coordinate, Z, for different values of the radiation number, N„ for the 

/v 

particular parameter values shown in the title of the figure. The oxidizer and fuel mass fractions 
at the respective walls are 7^=0. 6 and 7^=0. 8. In our subsequent analysis we keep the 
same set of oxidizer and fuel mass fractions, and vary the location, width and intensity of the 
radiative loss zone. The above set of ( 7^, Y FF ) represents a typical case and is employed 
extensively in the following analysis: the qualitative trends for other 7^ and Y FF values 
are similar. The thickness of the radiative loss zone is 0.04 for all values of N x and the 
separation distance of the loss zone from the stoichiometric flame location is zero. We observe 
that the flame temperature profile is uniformly lowered as the value of N R increases. Also, 
the flame temperature peak moves toward the fuel wall as the value of N R is increased. The 
drop in flame temperature, as well as the shift of the peak, become more prominent for higher 
values of N R . For a value of N R greater than 383, w-e do not obtain a steady state 
temperature profile, indicating the occurrence of a radiative extinction. This maximum, or upper 
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bound, for N R is then defined as A^ m>r ,. w i.e., for this case. We also 

note from Figure 4 that there is a change of slope of the temperature profile in the radiative loss 
zone for higher values of N^, i.e., between Z R _ and Z R 

Figure 5 shows the nondimensional reaction rate term ((l+<f>)Dr) for the same 
situation. We observe that the reaction rate profile collapses for increasing N R values. The 
reaction rate peak also moves towards the fuel side; this movement becomes more conspicuous 
for higher values of N R . We notice that the reaction rate profile has managed to move nearly 
beyond the rightmost side of the radiation loss zone (indicated by the dashed lines at Z R _ and Z R ) 
for the highest value of N R . 

We now focus on the temperature and species profiles for the situation when N R = 383 
for the above case, i.e., at the brink of extinction. Figure 6 also shows the temperature and 
species profiles for the same flame for an infinite reaction rate. We notice that when N R = 383, 
the slope of the y Q profile is quite different from its IRR counterpart. On the other hand, the 
slope of the y F profile follows the IRR y f profile closely until a Z-value of about 0.3, when 
its slope starts decreasing. This plot therefore demonstrates explicitly the contrast between the 
IRR situation and the finite chemistry situation with appreciable radiative losses. The migration 
of the peaks of temperature and reactivity profiles is striking. Also, an abrupt change of the 
temperature profile seems to take place in the zone of radiative losses, i.e. , between Z_ and Z R .. 
We add for emphasis that from the strictly physical viewpoint, the finite-rate solution has 
attained a rather extreme form, since the reaction zone has almost completely propagated through 
the loss zone. In Figure 6, we see that the loss zone is now on the oxidizer side of the reaction 
rate profile. As we shall see, extreme cases like this are not the norm. They are also physically 
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unrealistic, though mathematically permissible in our simplified model with a prescribed heat 
loss function. 

We illustrate the details of the flame structure in Figure 7, where we plot the 
contributions of the different terms in the energy equation when the steady state condition has 
been achieved. The loss term is given by (N / Js ( )sech 2 (B(Z-Z R )) and the diffusion term, as 
in equation (9), is (1 /s^x^. We have already noted from Figure 5 that for N R = 383 the 
reaction rate profile has penetrated through the radiative loss zone. Figure 7 indicates that the 
diffusion term recovers the radiative losses almost entirely and the reaction term doesn’t 
contribute to the diffusion term in such a recovery process. This represents a completely 
different physical problem, when the radiative loss term exists on the oxidizer side of the 
primary reaction zone (flame). This result is, as already mentioned, clearly in conflict with our 
hypothesis that the heat losses take place on the fuel side of the flame due to flame-generated 
particulates. This occurs because our hypothetical radiative loss profile is simply a prescribed 
function in Z, and as such, it does not contain any mechanisms for loss-zone movement as the 
temperature and species profiles change, as a real soot zone invariably must. 

In order to observe the effect of a thicker loss zone, we now increase A Z R to a value 
of 0.1, see Fig. 8. We notice that the drop in the temperature profile is more significant in this 
case and the flame extinguishes at a lower value of the radiation number, viz., for N^=132. 
We use this opportunity to note that (A/^AZ^ is approximately 14.2 for the first case and 
13.2 for the second case. As shown in Wichman and Ray (14) for the simple top-hat loss 
profile, it appears the extinction results are best correlated with the function N R &Z R , although 
the proper method of evaluating A Z R is not as straightforward as our estimate suggests. The 
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correct method of evaluating A Z R is discussed later in the article. 

Next, we consider the situation when the left-most side of the loss zone is sufficiently 
removed from Z f for a flame with >' oo =0.6 and ^=0.8. The thickness of the loss zone 
is AZ^=0.06 and the separation distance, A, is 0.1 in this case. Figure 9 shows that the 
flame temperature decreases with increasing N R . Here, the movement of the peak 
nondimensional flame temperature is not pronounced, though it does move toward the fuel side. 
Correspondingly, Figure 10 shows the variation of the reaction term, (1 +<$>)Dr, for increasing 
values of Nj,. As mentioned for the preceding cases, therefore, the reaction zone does not 
always propagate through the loss zone. A sufficient separation and magnitude of the loss term 
appear sufficient to block the through-transit. 

We note that the reaction rate peak is always to the left of the temperature peak, i . e. , 
Z f <Z<Z„. This is in accordance with the results obtained for pure diffusion flames without 
radiative loses [9], as discussed before. 

Figure 1 1 is an extinction plot for the case w'hen AZ^=0.06 and A=0. Extinction 
values of N R are plotted as a function of Z ^ the theoretical flame location in the mixture 
fraction coordinate. We notice that for a given value of the oxidizer mass fraction at the wall, 
(^gXr mjutwn increases as Z^ is decreased. A decrease in Zj implies an increase in 4>, 
which, for a given Y^, produces an increase in Y fr As Y FF increases, the reaction rate 
becomes more vigorous and it becomes more difficult to extinguish the fiame through radiative 
losses. This explains the nature of the curves that we obtain on the extinction plot. Also, for 
the same value of Z„ i.e., for the same value of <{>, a lower value of Y ^ indicates a 
correspondingly smaller value of Y Fr and hence, the reaction rate also becomes smaller in 
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magnitude. It is then easier to extinguish the flame. This explains why the curves in Figure 1 1 
all shift toward the left for decreasing values of Y^. 

Our focus is next shifted to some quantities of practical interest. We evaluate the heat 
transfer to the wall from flames with the same stoichiometry (T 00 = 0.6 and Y FF =0. 8) but 
with different thicknesses of the radiative loss zones and for different separation distances A 
from Z f Let Q wo denote the heat transferred to the oxidizer wall by the flame per unit 
surface area of the wall. We reckon that the oxidizer wall will have a stronger effect on the 
flame than the fuel wall owing to the proximity of the flame to the oxidizer wall. The flame 
transfers heat to the oxidizer wall by means of both conduction and radiation, and hence, 
Qw.o = Qw.o.cor^Qw.o /a s’ where the conduction flux is Q w _ 0 .co^~ x ^ dT l dx )\ x .L the 

radiative flux is Q WO/aJ =0.5'/.^\dq F jds)dx. We assume that half of the radiative losses travel 
to each wall: this assumption is reasonable in the thermally-thin limit we consider here. We 
can transform the expressions for Q WOcond and Q VOtrad to the Z coordinate and normalize Q wo 
by the reference conductive flux X o iT f -T^IL. The normalized 
I : .o + 0- 5x (1/^^/q (l/^o)^^ 2 )^- The quantity Q wO s 0 is plotted in 
Figure 12. It is apparent that the heat transfer characteristics do not depend strongly on the 
separation distance A, and consequently, we see four reasonably distinct groups of curves 
corresponding to loss zones of four different thicknesses. However, as is evident from the plot, 
the separation distance A does become important for higher values of close to 

extinction. We will notice that, nearing extinction, the flame attempts to reduce the heat losses 
to the wall as much as possible. Also, the value of N R required for extinction is higher when 
the heat loss zone is very thin, as intuitively obvious. 
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We also plotted Q vo s 0 as a function of the quantity N R (2/B). Figure 13 clearly 

shows that the quantity N R (2/B), which is approximately the value of the integral 

f^Njech^BCZ-Z^dZ, is able to collapse the wall heat transfer data except very near 

extinction. Thus, when plotted against N R (2jB), Q W( ^ 0 does not reveal any appreciable 

dependence on either the separation distance A or even the thickness of the loss zone A Z R . 

The correlation in the straight-line region is given approximately by Q w or s Q “\ .9{N R i^ll B)) .2 , 

which is reasonably accurate until the curves make their Final turn toward the abscissa. 

Another quantity of practical interest is the radiative fraction y, given by the ratio 

The quantity q ^ is the integral of the radiative loss term, i.e., 

J o N^ech 2 [B(Z-Z R )]dZ and q Total is the integrated value of the reaction rate in mixture 

fraction space, i.e., J (1 +<$>)DrdZ. From Figure 14 we notice that Q TotcJ decreases with 

increasing values of N R . This happens because with increased intensity of the radiative loss 

zone, reaction rate values decrease as already observed in Figures 5 and 10. For thicker loss 

zones, the drop in q ToccJ with increasing values of N R is more rapid. 

We have already noted in section 2.2 that the integral of the radiative loss term profile 

is approximately (N^s^il/B). Hence, it is of interest to plot the total heat release q Jgtc! as 

a function of the quantity N R (2/B). The result is shown in the Figure 15. Figure 15 indicates 

•• 

that the quantity N R (2/B) characterizes the total heat release rate very well and the curves for 
different loss zone thicknesses virtually collapse on one another except for large values of N R 
close to extinction. The correlation in the linear region is given by q TouzJ =-\5.25N R (2/B) + \6&. 
The constants in this formula depend on the global stoichiometry. 

Figure 16 illustrates the variation of % as a function of N R for different thicknesses 
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of the loss zones and for A =0. We observe that x increases with increasing N R for a 
flame with a given loss zone thickness. The integrated quantity q ^ increases with N R 
and, since correspondingly the q Total values decrease, x. which is a ratio of the above 
quantities, increases. In order to produce a given value of x. a higher value of N R is 

required for a flame with a thinner loss zone. Similar to the study of q Toul p we plot x as 

a function of N R (2/B) in Figure 17. It is clear from the Figure that the use of N R (2/B ) 
collapses the data very well except close to extinction. Here, the correlation in the linear 
segment is x = 0-038iV A (2/3); once again, the multiplicative constant must be a function of 
global stoichiometry. Figure 18 shows the variation of the drop in flame peak temperature, 
At^, as a function of the radiative fraction x- If we denote the maximum temperature by 
then Axy is defined as 1-t^. We recall that the temperature has been normalized in such 

a way that the peak nondimensional temperature for the infinite reaction rate situation always 

has the value of unity, regardless of the oxidizer and fuel mass fractions. Thus, At. 
represents the drop in peak temperature for finite rate chemistry and radiative loss situation, in 
comparison to the IRR situation. The increase in At^ with x was almost linear for smaller 

values of x. with At^x+O l as the correlating function. However the curves for the 

different loss zone thicknesses diverged from one another for higher values of x- 

It has been previously mentioned (section 2.2) that the results for the sechr heat loss 

profile can be compared with the results of Wichman and Ray [14] for the top-hat profile. 

However, the thickness A Z R of the top-hat profile must be chosen to be 2/3, where the value 

of 3 is determined from the choice of the thickness of the sec } r profile. For example, when 

is chosen to be 0.06, the constant 3=99.7 and consequently 
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_ te =2/99.7“0.02. As shown in [14], both analytical and numerical methods were used 
to determine the extinction value of N R for the top-hat profile. Here, we compare analytical 
and numerical results for the top-hat profile with the numerical solutions for the sech : profile. 
Figure 19 depicts the extinction N R values plotted as a function of Z { when F oo =0. 7, 
(AZ^ cA i= 0.06, (AZ /t ) £op _ /u2/ =0.02 and A =0.1. The direction of increasing Y FF has also 

been indicated on the plot. The numerical solutions reveal that the sech 2 and the top-hat profiles 
produce very similar values. This indicates that the integrated value of the radiative 

loss term is the quantity which determines the extinction N R value. The extinction M R 
values obtained by analytical method are quite different from the numerical solution. However, 
inspection of the curves depicted in Figure 19 shows that the ratio of the analytically obtained 
values to the numerical solution is very nearly (±1.5%) 3.8 for all the Z f values plotted in 
Figure 19, exactly as in [14], This indicates that a simple modification of the analytical formula 
of [14] based on the inclusion of a correction factor should yield close correspondence between 
the analytical and numerical results. Thus, we use 




C 


Zj 1-0/2) p 




4D 0 Z/(1-Z / ) 3 


(M ) 3 


(17) 


where D o is the Damkohler number defined under Equation (1 1), b? is the extinction value 
of the reduced Damkohler number (see [14]) and 0=Z^+Z J? _, which we write as 0-2 Z r 
after using Z Rt =Z R +AZ^2, Z R _=Z R -&Z R /2. In other words, Z R is the value of the heat loss 
zone peak. With C= 1/3.8 this formula correlates the numerically-derived data of Fig. 19 within 
line width. This multiplicative factor will depend on the global stoichiometry, of course. 
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5. 


CONCLUSIONS 


We investigated the influence of a simple and hypothetical heat loss zone on a pure 
diffusion flame in detail. The loss profile was of the form of a sechr given by equation (12) and 
we varied the intensity and the width of the loss zone to study the diffusion flame response. The 
loss zone was postulated to lie always on the fuel side of the ideal Burke-Schumann flame. The 
location of the loss zone on the fuel side relative to the ideal flame location was also 
varied. In all situations, the increase of the radiation number N R results in a movement of 
the flame toward the fuel side. We found that for thin loss zones located close to Z ^ the 
reaction zone may even migrate through the flame to the fuel side of loss zone for significantly 
high values of N R . In such a situation the loss zone now lies on the oxidizer side of the flame, 
contrary to our initial postulate. This happens because our hypothetical loss zone is static and 
contains no mechanism for movement. The reaction rate profile, on the other hand, is free to 
move and hence locates itself on the fuel side of the loss zone in certain cases. 

Extinction plots were generated for different flames for given loss zone thicknesses 
(AZ^) and given separation distances (A). The plots indicated that for a given Y oa , an 
increase in Y FF results in higher values of extinction radiation number, N R . 

Nondimensional heat transfer rates to the oxidizer wall were also investigated. The 
results indicated that the separation distance A did not have a significant influence on the wall 
heat transfer characteristics. 

The total heat release in the combustion process, q To[al , was found to decrease with 
increasing values of N R and the rate of decrease was quite rapid for thicker loss zones. It was 
found that the quantity N R (2/B) characterizes q Total very well and the curves for the 
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different loss zone thicknesses and separation distances all collapse onto one another except near 
extinction. A correlation was produced in the linear region. We expect that the correlation 
constants depend upon the flame stoichiometry. 

Investigation of the radiative fraction showed that x increases with increasing values 
of N r and the rate of increase is steeper for thicker loss zones. The flame, however, 
extinguishes at a larger value of x for loss zones which are relatively thin. The quantity 
N r (2/B) collapses the q 7otcl and radiative fraction values very well except near flame 
extinction. A linear correlation formula was produced whose coefficients vary with global flame 
stoichiometry. The detailed determination of these coefficients will be the subject of a future 
work. The decrease in the flame temperature was nearly linear with radiative loss fraction, as 
shown in Fig. 18. The linearity, away from near-extinction, of all of these results is very 
encouraging for the development of a simplified description of soot radiation in flames. Once 
a suitable N R is defined and an estimate is made of the soot layer thickness 2 IB in mixture- 
fraction space, correlations resembling those of Figs. 13, 15, 17, 18 can be generated for the 
important overall heat-transfer quantities. The estimation of the soot-layer thickness will be the 
subject of a future work. 

Finally, we demonstrated that the extinction formula derived in Ref. [14] could be easily 
modified to suit this sechr model by replacing the top-hat loss zone thickness with 2/5 for large 
5. The hop-hat and sech 2 profiles thus modified produce exactly the same N KeBijiaion vs 
curve, see Fig. 19. The correlation with the analytical formula derived in [14] is given by 
equation (17) with the new effective loss-zone thickness 2/5. We expect that the multiplicative 
empirical constant (here 3.8), will depend fairly weakly on the global stoichiometry. Many 
cases must be examined in order to determine this dependence. 
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NOMENCLATURE 


a Constant in asymptotic formulation defined by Eq. (35. i) 

a Planck mean absorption coefficient 

A Pre-exponential factor 

b Reduced Damkohler number defined by Eq. (35. iii) 

b' Reduced Damkohler number with zero heat losses 

c 0(1) constant in asymptotic formulation, see (Eq. (35. ii) 

DJD o Damkohler number; Damkohler number with no heat losses 
D gr D 1 e Diffusion coefficients for oxidizer, fuel 

E Activation energy 

f v Soot volume fraction 

E Specific enthalpy of species / 

h Rescaled enthalpy loss, h = $\H(Z R _)\ 

H Enthalpy defect, H=x +y o *y F ~\ 

H o Enthalpy defect with zero heat losses, // =0 

L Combined heat loss; distance between fuel and oxidizer walls 

Radiation number, see Sec. 2.4 

q R Nondimensional radiant energy flux, see (Eq. (1) and above Eq. (2) 
Q F Heat release per unit mass of fuel 

Q f Nondimensional heat release 

t Nondimensional reaction term, -r =y & y^xp[ — P(1 — x)/[l — cc(l - t)]) 

R Universal gas constant 

s 7 s g Mass coordinate, s= f X pdx, s g = f L pdx 

J O Jo 



s Nondimensional mass coordinate, s =s/s 

S Redefined dependent variable, S=x-H-l-(y g +y F ) 

T,TpT g ,T R Temperature (flame, ambient, radiation zone) 

u Flow velocity 

U g (Z) Heaviside step function 

w Reaction term, w=pAY g Y f £xp(-E/RT) 

x Spatiai coordinate 

y F ,y o Rescaled fuel, oxidizer mass fractions 

Z^Y g Fuel, oxidizer mass fractions 

Z Physical coordinate, Z=l-s, coincides with mixture fraction 

Z R ^YZ R JL R Z-values at right, left, middle of radiation-loss zone 


GREEK 


a Enthalpy ratio, a =1 -7 JT, 

P Zeldovic'n number, P=a£/R7^ 

A Z B Radiation-loss zone thickness, AZ =Z, -Z. 
e Small dimensionless parameter 

t| Rescaled mixture fraction, see above Eq. (31) 

0 Sum of Z R „ and Z R _, Q=Z R ^Z R _ 

X Thermal conductivity 

; Nondimensional physical coordinate, i=s 

p Density 
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o Stefan-Boltzmann constant 

t NondimensionaJ temperature, 

4> Global equivalence ratio 




$ 


Scaled value of 5 in reaction zone, see above Eq. (31) 
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LIST OF CAPTIONS FOR FIGURES 


Figure 1: 

Figure 2: 

Figure 3: 
Figure 4: 

Figure 5: 
Figure 6: 

Figure 7: 
Figure 8: 

Figure 9: 

Figure 10: 

Figure 11: 

Figure 12: 
Figure 13: 


The problem geometry, including the diffusion flame (DF), the radiation (soot) 
zone and the porous diffusive walls at x = 0 and x=L. The reactant influx is 
purely diffusive. 

Flame temperatures via simple correlation. Note that the maximum for 

Y FF =0.25, 0.3 occurs near 7^=0. 4: we do not use these curves for this very 

reason. 

The sech 2 heat-loss profile showing the separation A between Zj and Z ? _, 
and our (preliminary) definition of A Z R (which we shall subsequently modify). 

Influence of increasing N R on the temperature distribution in a diffusion flame. 
Here the flame penetrates the loss zone, which is unrealistic unless oscillations 
occur. When N R >3$3, r crashes everywhere to zero. 

Same as Fig. 4 for the reaction rate. Here the migration of the flame is 
immediately obvious. 

Profiles of t, y ^ y F and reaction rate for infinite and finite-rate cases, the 
latter just prior to extinction. Note that the y F profile is virtually unchanged 
from its IRR value until inside the reaction zone. 

The flame structure for Fig. 4 at the brink of extinction. 

Influence on temperature field of increasing N R on the same DF as Fig. 4 
except the loss zone is 2.5x thicker, AZ /J =0.1. 

Influence of non-zero initial displacement A between the loss zone and the DF 
(at Z1 Here, the flame does not penetrate the loss zone before extinction. 

Same as Fig. 9 for the reaction rate. The right-ward movement in mixture 
fraction space is obvious. 

Extinction plot for AZ /? =0.06 and zero initial displacement, A=0. Along 
each curve, Y FF increases as Z { decreases. The largest Y FF values have 

the highest (A^)^ values. 

Total heat flux (nondimensional) to the oxidizer wall versus N R . 

A replot of Fig. 12 with new abscissa N R (2/B). The correlation of the straight- 
line region is (heat flux) = 1.9/V ff (2/fi) + 5.3. 



Figure 14: The total heat flux as a function of N R for various loss-zone thicknesses and 

zero displacement, A =0. As N x increases, the integrated reaction rate 
decreases monotonically. 

Figure 15: Replot of Fig. 14 showing the collapse of the data to a single line except very 

near extinction. 

Figure 16: Radiant flux fraction versus N R . 

Figure 17: A replot of Fig. 16 with abscissa N R (2/B) showing collapse of data to a single 

curve. 

Figure 18: The decrease in flame temperature versus radiant flux fraction shows a linear 

functional form that is virtually independent of the loss-zone width, AZ^,. 

Figure 19: Plot of versus Z / for the seek 2 loss profile and the previously-examined 

top-hat loss profile [14], The close agreement between the numerical solutions 
indicates that the correlation with the analytical formula will be outstanding when 
the latter formula is multiplied by (3.8)'* = 0.26 1 . 



TABLE 1: FLAME TEMPERATURE VERSUS DATA 
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